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Abstract. We propose the graph description of Teichmiiller theory of surfaces with marked 
points on boundary components (bordered surfaces). Introducing new parameters, we formu- 
late this theory in terms of hyperbolic geometry. We can then describe both classical and 
' quantum theories having the proper number of Thurston variables (foliation-shear coordi- 

nates), mapping-class group invariance (both classical and quantum), Poisson and quantum 
algebra of geodesic functions, and classical and quantum braid-group relations. These new 
algebras can be defined on the double of the corresponding graph related (in a novel way) 
to a double of the Riemann surface (which is a Riemann surface with holes, not a smooth 
Riemann surface). We enlarge the mapping class group allowing transformations relating 
different Teichmiiller spaces of bordered surfaces of the same genus, same number of boun- 
^-j- , dary components, and same total number of marked points but with arbitrary distributions 

;> of marked points among the boundary components. We describe the classical and quantum 

s ! . algebras and braid group relations for particular sets of geodesic functions corresponding 

to A n and D n algebras and discuss briefly the relation to the Thurston theory. 

C> Key words: graph description of Teichmiiller spaces; hyperbolic geometry; algebra of geo- 

T ~' ■ desic functions 
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1 Introduction 

Recent advances in the quantitative description of the Teichmiiller spaces of hyperbolic structures 
were mainly based on the graph (combinatorial) description of the corresponding spaces [191 [7]. 
The corresponding structures not only provided a convenient coordinatization together with the 
mapping class group action, they proved to be especially useful when describing sets of geodesic 
functions and the related Poisson and quantum structures [3]. Combined with Thurston's theory 
of measured foliations [22[ [20] , it led eventually to the formulation of the quantum Thurston 
theory [5]. The whole consideration was concerning Riemann surfaces with holes. A natural 
generalization of this pattern consists in adding marked points on the boundary components. 
First, Kaufmann and Penner [15] showed that the related Thurston theory of measured foli- 
ations provides a nice combinatorial description of open/closed string diagrammatic. Second, 
if approaching these systems from the algebraic viewpoint, one can associate a cluster algebra 
(originated in [10] and applied to bordered surfaces in [11]) to such a geometrical pattern. 

The aim of this paper is to provide a shear-coordinate description of Teichmiiller spaces of 
bordered Riemann surfaces, to construct the corresponding geodesic functions (cluster variables), 
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and to investigate the Poisson and quantum relations satisfied by these functions in classical 
case or by the correspondent Hermitian operators in the quantum case. 

In Section [21 we give a (presumably new) description of the Teichmuller space of bordered (or 
windowed) surfaces in the hyperbolic geometry pattern using the graph technique supporting it 
by considering a simplest example of annulus with one marked point. It turns out that adding 
each new window (a new marked point) increases the number of parameters by two resulting 
in adding a new inversion relation to the set of the Fuchsian group generators. We explicitly 
formulate rules by which we can construct geodesic functions (corresponding to components of 
a multicurve) using these coordinates; the only restriction we impose and keep throughout the 
paper is the evenness condition: an even number of multicurve lines must terminate at each 
window. 

In Section [3l we construct algebras of geodesic functions postulating the Poisson relations 
on the level of the (old and new) shear coordinates of the Teichmuller space. We construct 
flip morphisms and the corresponding mapping class group transformations and find that in 
the bordered surfaces case we can enlarge this group allowing transformations that permute 
marked points on one of the boundary components or transfer marked points from one compo- 
nent to another thus establishing isomorphisms between all the Teichmuller spaces of surfaces of 
the same genus, same number of boundary components, and the same total number of marked 
points. In the same section, we describe geodesic algebras corresponding (in the cluster termi- 
nology, see [11]) to A n and D n systems. Whereas the A n -algebras have been known previously 
as algebras of geodesies on Riemann surfaces of higher genus [161 [P7] (their graph description in 
the case of higher-genus surfaces with one or two holes see in [1]) or as the algebra of Stockes 
parameters [6j [23] , or as the algebra of upper-triangular matrices [2] , the D n -algebras seem to 
be of a new sort. Using the new type of the mapping class group transformations, we prove the 
braid group relations for all these algebras. 

Section [4] is devoted to quantization. We begin with a brief accounting of the quantization 
procedure from [3] coming then to the quantum geodesic operators and to the corresponding 
quantum algebras. Here, again, the quantum D n -algebras seem to be of a new sort, and we 
prove the Jacobi identities for them in the abstract setting without appealing to geometry. We 
also construct the quantum braid group action in this section. 

In Section [5l we describe multicurves and related foliations for bordered surfaces, that is, we 
construct elements of Thurston's theory. There we also explicitly construct the relevant doubled 
Riemann surface, which, contrary to what one could expect, is itself a Riemann surface with 
holes (but without windows). We transfer the notion of multicurves to this doubled surface. 
Note, however, that the new mapping class group transformations, while preserving the mul- 
ticurve structure on the original bordered surface, change the topological type of the doubled 
Riemann surface, which can therefore be treated only as an auxiliary, not basic, element of 
the construction. Using this double, we can nevertheless formulate the basic statement similar 
to that in [5], that is, that in order to obtain a self-consistent theory that is continuous at 
Thurston's boundary, we must set into the correspondence to a multicurve the sum of lengths of 
its constituting geodesies (the sum of proper length operators in the quantum case). In the same 
Section [5l we describe elements of Thurston's theory of measured foliations for bordered Rie- 
mann surfaces and the foliation-shear coordinate changings under the "old" and "new" mapping 
class group transformations. 

We discuss some perspectives of the proposed theory in the concluding section. 

We tried to make the presentation as explicit and attainable as possible, so there are many 
figures in the text. 
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2 Graph description and hyperbolic geometry 

2.1 Hyperbolic geometry and inversion relation 

2.1.1 Graph description for nonbordered Riemann surfaces 

Recall the graph description of the Fuchsian group, or the fundamental group of the surface S S)S) 
which is a discrete finitely generated subgroup of PSL(2,M). We consider a spine T 9tS corre- 
sponding to the Riemann surface £ 9jS with g handles and s boundary components (holes). The 
spine, or fatgraph T giS is a connected graph that can be drawn without self-intersections on 
has all vertices of valence three, has a prescribed cyclic ordering of labeled edges entering each 
vertex, and is a maximum graph in the sense that after cutting along all its edges, the Riemann 
surface decomposes into the set of polygons (faces) such that each polygon contains exactly one 
hole (and becomes simply connected after plumbing this hole). Since a graph must have at least 
one face, we can therefore describe only Riemann surfaces with at least one hole, s > 0. The 
hyperbolicity condition also implies 2g — 2 + s > 0. We do not impose restrictions, for instance, 
we allow edges to start and terminate at the same vertex, allow two vertices to be connected 
with more than one edge, etc. We however demand a spine to be a cell complex, that is, we do 
not allow loops without vertices. 

Then, we can establish a 1-1 correspondence between elements of the Fuchsian group and 
closed paths in the spine starting and terminating at the same directed edge. Since the terms 
in the matrix product depend on the turns in vertices (see below), it is not enough to fix just a 
starting vertex. To construct an element of the Fuchsian group A ff . s , we select a directed edge 
(one and the same for all the elements; see the example in Fig. H] where it is indicated by a short 
fat arrow) , then move along edges and turns of the graph without backtracking and eventually 
turn back to the selected directed edg^E 

We associate with the ath edge of the graph the real Z a and set [7] the matrix of the Mobius 
transformation 

-e Za/2 



X z a = [ „- Za/ 2 (2-1) 



each time the path homeomorphic to a geodesic 7 passes through the ath edge. 

We also introduce the "right" and "left" turn matrices to be set in the proper place when 
a path makes the corresponding turn, 

R = ( -1 ) ' L = r2 = ( -1 -1 

and define the related operators Rz and L^, 



R z = RX Z 
Lz = LXz - 



e -Z/2 _ e Z/2 

e z / 2 

_ e -Z/2 e Z/2 



An element of a Fuchsian group has then the structure 
P 7 = LX Zn RX Zn _ r • • • RX Z2 RX Zl , 



1 If the last edge was not the selected one but its neighboring edge, the very last move is turning to the selected 
edge, that is, we add either R- or L-matrix; if the last edge coincides with the selected edge, we do not make the 
last turn through the angle 2tt. This results in the ambiguity by multiplication by the matrix —Id = R A , but it 
is inessential as we deal with the projective transformations. 
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and the corresponding geodesic function 

G 7 = trP 7 = 2cosh(^ 7 /2) (2.2) 
is related to the actual length £ 7 of the closed geodesic on the Riemann surface. 

2.1.2 Generalization to the bordered surfaces case 

We now introduce a new object, the marking pertaining to boundary components. Namely, we 
assume that we have not just boundary components but allow some of them to carry a finite 
number (possibly zero) of marked points. We let <5j, i = 1, . . . , s, denote the corresponding num- 
ber of marked points for the ith boundary component. Geometrically, we assume these points 
to lie on the absolute, that is, instead of associating a closed geodesic to the boundary compo- 
nent in nonmarked case, we associate to an ith boundary component a collection comprising 5i 
infinite geodesic curves connecting neighbor (in the sense of the surface orientation) marked 
points on the absolute (can be the same point if 5i = 1) in the case where 5i > 0. All these 
additional geodesic curves are disjoint with each other and disjoint with any closed geodesic on 
the Riemann surface. In [15], these curves were called windows. We denote the corresponding 
windowed surface S s < 5, where 

S = {S 1 ,...,5 S } (2.3) 

is the multiindex counting marked points on the boundary components {5i can be zero) whereas 
s is the number of boundary components. We call such Riemann surfaces the windowed, or 
bordered Riemann surfaces. 

Restrictions on g, s, and the number of marked points #<5 can be uniformly written as s > 



#s+i 



> 0, that is, we allow two new cases g = 0, s = 1, #5 > 3 and g = 0, 



and 2g - 2 + s + 
s = 2, #5> 1. 

We want now to generalize the graph setting to the case where we have boundary components 
with marked points. However, as the example below shows, in order to define inambiguously the 
corresponding hyperbolic geometry, when introducing a marked point on the boundary, we must 
simultaneously introduce one more additional parameter. This is because, as we shall demon- 
strate, introducing a marked point adds a new inversion relation that preserves the orientation 
but not the surface itself, that is, we invert a part of the Riemann surface through a boundary 
curve without taking care on what happen to the (remaining) part of the surface because, in our 
description, this (outer) part is irrelevant. Such an inversion relation leaves invariant the new 
added geodesic, that is, the window. However, there is a one-parameter family of such inversions 
for every window, and in order to fix the ambiguity we must indicate explicitly which point on 
the new geodesic is stable w.r.t. such an inversion. Recall that because of orientation preser- 
vation, two ends (on the absolute) of this new geodesic must be interchanged by the inversion 
relation; it is therefore a unique point that is stable. 

We describe this situation by considering new types of graphs with pending vertices. Assume 
that we have a part of graph having the structure as in Fig. [TJ 

Then, if a geodesic line comes to a pending vertex, it undergoes the inversion, which stems 
to that we insert the inversion matrix F, 



1 
-1 



(2.4) 



into the corresponding string of 2 x 2-matrices. For example, a part of geodesic function in Fig. [T] 
that is inverted reads 



• Xy 1 LXzFXzLXy 2 ■ ■ ■ > 
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Figure 1. Part of the spine with the pending vertex. The variable Z corresponds to the respective 
pending edge. Two types of geodesic lines are shown in the figure: one that does not come to the edge 
Z is parameterized in the standard way, the other undergoes the inversion with the matrix F (|2.4p . 

whereas the other geodesic that does not go to the pending vertex reads merely 

• • • Xy 1 RXy 2 ■ ■ ■ ■ 

We call this new relation the inversion relation, and the inversion element is itself an element of 
-PSX(2,]R). We also call the edge terminating at a pending vertex the pending edge. 
Note the simple relational, 

XzFXz = X2z- 

We therefore preserve the notion of the geodesic function for curves with inversions as well. 
We consider all possible paths in the spine (graph) that are closed and may experience an 
arbitrary number of inversions at pending vertices of the graph. As above, we associate with 
such paths the geodesic functions (here, we let Zi denote the variables of pending edges and Yj 
all other variables) 

G 7 = tr P 7 = 2 cosh(V2) = tr LX Zn FX Zn RX Yn _^ ■ ■ ■ RX Zl FX Zl . (2.5) 
We have that, for the windowed surface £ 9i <5, the number of the shear coordinates Z a is 

s 

#Z Q = 6< 7 -6 + 3s + 2^5 i , 

i=i 

and adding a new window increases this number by two. 

Before describing the general structure of algebras of geodesic functions, let us clarify the 
geometric origin of our construction in the simplest possible example. 

2.2 Annulus with one marked point 

The simplest example is the annulus with one marked point on one of the boundary components 
(another example of disc with three marked points will be considered later). Here, the geometry 
is as in Fig. [2] where we let the closed line around the neck (the blue line) denote a unique 
closed geodesic corresponding to the element Pi of the Fuchsian group to be defined below, the 
winding to it line (the red line) is the boundary geodesies from the (ideal) triangle description 
due to Penner and Fock, and the lower geodesic (the green line) is the new line of inversion (the 
window). We indicate by bullet the stable point and by cross the point of the inversion line that 
is closest to the closed geodesic. 

2 In particular, we would consider inversion generated by a Mobius element type, (|2.ip . say, F = Xw not just 
F — Xo- But then XzXwXz — X 2 z+w, so we can always adsorb W into Z thus producing no new factors; we 
therefore stay with our choice of F . 
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Figure 2. Geodesic lines on the hyperboloid: dotted vertical line is the asymptote going to the marked 
point on the absolute, closed blue line is a unique closed geodesic; red line is the line from the ideal 
triangular decomposition asymptotically approaching the asymptote by one end and the closed geodesic 
by the other; green line is the line of inversion whose both ends approach the marked point; we let the 
bullet on this line denote the unique stable point under the inversion and the cross denote the point that 
is closest to the closed geodesic. 



The same picture in the Poincare upper half-plane is presented in Fig. [3l There, the whole 
domain in Fig. [2] bounded below by the bordered (green) geodesic line and above by the neck 
geodesic (blue) line is obtained from a single ideal triangle with the vertices {e z+Y , oo, 0} upon 
gluing together two (red) sides of this ideal triangle. We now construct two (hyperbolic) ele- 
ments: Pi that is the generating element for the original hyperbolic geometry and the new 
element P\\ that corresponds to the inversion w.r.t. the lower (green) geodesic line in Figs. [2] 
and [3l Adding this new element obviously changes the pattern, but because the Fuchsian pro- 
perty retains, the quotient of the Poincare upper half-plane under the action of this new Fuchsian 
group must be again a Riemann surface with holes. As we demonstrate below, this new Riemann 
surface is just the double of the initial bordered Riemann surface. 

For this, we use the graph representation. The corresponding fat graph is depicted in Fig. [H 
This graph with one pending edge and another edge that starts and terminates at the same 
vertex is dual to an ideal triangle {e z+Y , oo, 0} in which two (red) sides are glued one to another 
(the resulting curve is dual to the loop) and the remaining (green) side is the boundary curve 
(dual to the pending edge). We mark the starting direction by the fat arrow, so the element Pi is 

/ e -Y/2 i e Y/2 _ Z+Y/2 \ 

Pi = XzLXyLXz = I e _ z Z Y/2 Q . (2.6) 

Apparently, the corresponding geodesic function G\ is just e~ Y / 2 + e y / 2 , so the length of the 
closed geodesic is \Y\ as expected. 

We now construct the element Pn- Note that this element makes the inversion w.r.t. the 
geodesic between and oo, so we set the matrix F first (since the multiplication is from right 
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Figure 3. The hyperbolic picture corresponding to the pattern in Fig.O preimages of red boundary line 
are red half-circles, preimages of the inversion line are green half circles (the selected one connects the 
points oo and on the absolute), and the preimage of the closed geodesic is the (unique) blue half-circle; 
the points e z and e z+Y on the absolute are stable under the action of the corresponding Fuchsian element 
Pi (|2.6p ; the bullet symbols are preimages of the point that is stable upon inversion (the one that lies 
on the geodesic line between oo and is i in the standard coordinates on the upper half-plane) and the 
dotted half-circles connect the point e z+Y with its images (one of which is — e~ z ~ Y ) under the action 
of the inversion element F. We also mark by cross the point ie z+Y < 2 of the green geodesic line that is 
closest to the closed geodesic. The invariant axis of the new element P\\ (|2.7[) and some of its images 
under the action of (|2.6|) are depicted as cyan half-circles; £2 is from (|2.8|) . 



to left, this matrix will be rightmost). Then, the rest is just the above element P\\ 

( e z+Y ' 2 \ 

Pn = X Z LX Y LX Z F = P,F = ^ e _ Y/2 + &Y/2 g _ z _ y/2 J , (2.7) 

and the corresponding geodesic function G\\ is 2 cosh(Z+Y/2) so the length of the corresponding 
geodesic (but in a geometry still to be defined!) is \2Z + Y\. 

We now consider the action of these two elements in the geometry of the Poincare upper half- 
plane in Fig. [3j It is easy to see that the element Pi has two stable points: e z (attractive) and 
e z+Y (repulsive). It also maps 00 — > 0, e z + e z+Y — > 00, etc. thus producing the infinite set of 
preimages of the red geodesic line in Fig. [2] upon identification under the action of this element. 

The element F first interchanges and 00 and e z+Y and — e~ z ~ Y thus establishing the 
inversion (inversion) w.r.t. the green geodesic line. The only stable point of this inversion is the 
point of intersection of the two above geodesic lines, and it is the point i in the upper complex 
half-plane for every Z + Y '. Further action is given by P\ and, in particular, it maps 00 back 
to 0, so £1 = is a stable point of Pn. Another stable point is 

e Z+Y/2 _ e -Z-Y/2 

^ = e Y/2 + e -y/2 ' ( 2 ' 8 ) 

and it is easy to see that the two invariant axes of Pi and Pn never intersect. Adding the 
element Pn to the set of generators of the new extended Fuchsian group we therefore obtain 
a new geometry. 

First, let us consider the special case where the stable point on the inversion curve coincides 
with the point that is closest to the closed geodesic. Then, apparently, the inversion process 
exhibits a symmetry depicted in Fig. [5j Considering the Riemann surface depicted in Fig. [2 
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Figure 4. The graph for annulus with one marked point on one of the boundary components. Examples 
of closed geodesies without inversion (I) and with inversion (II) are presented. The short fat arrow 
indicates the starting direction for elements of the Fuchsian group. 




Figure 5. The doubled Ricmann surface obtained upon inversion w.r.t. the green geodesic in the case 
where the stable point coincides with the point closest to the closed geodesies (blue line) (the cross then 
coincides with the bullet). 

we chop out all its part that is below the green (inversion) line. We then obtain the double of 
the Riemann surface merely by inverting it w.r.t. the green line taking into account the obvious 
(mirror) symmetry that takes place in this case. We then obtain from the hyperboloid with 
marked point at the boundary component the sphere with two identical cycles (images of the 
closed geodesic) and one additional puncture (hole of zero length) , as shown in Fig. [5j 

What happens if, instead of the stable point marked by cross, we have arbitrary stable point 
(bullet in Figs. [2] and [3])? Actually, we can answer this question just from the geometrical 
standpoint. Indeed, since in the pattern in Fig. El points on the inversion geodesies that lie 
to both sides from the asymptote are close, they must remain close in the new geometry. But 
the image of each such point is shifted by a distance that is twice the distance D (along the 
inversion line, which is a geodesic line) between the stable point (bullet) and the symmetric 
point (cross). This means that, in the new geometry, the points on the inversion line separated 
by a distance 2D must be asymptotically close as approaching the absolute in the pattern of 
Fig. [3l This means in turn that the corresponding geodesic in the new geometry is just a geodesic 
approaching the new closed geodesic of length 2D. 

It remains just to note that, from the pattern in Fig. [3l 

D = \Z + Y/2\, 
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Figure 6. The doubled Riemann surface obtained upon inversion w.r.t. the green geodesic in the case 
where the stable point (marked by •) is arbitrary. The closed in the asymptotic geodesic sense points in 
the new geometry are those on different coils of the spiraling green geodesies, which has the asymptotic 
form of the double helix. The separation length is asymptotically equal to \2Z + Y\. The cyan line is the 
new closed geodesic (the invariant axis of the element Pn). We let two crosses denote the points on the 
inversion line that are closest to the two copies of the initial closed geodesic line; the geodesic distance 
between them is also \2Z + Y\. 

that is, the perimeter of the new hole is \2Z + Y\, and it coincides with the length of the 
new element P\\ (|2.7p . which is therefore the element of the new, extended, Fuchsian group 
corresponding to going round the new hole. In Fig. [6j we depict this new geometry. It is also 
interesting to note that we now again, as in the symmetrical case, have two (homeomorphic) 
images of the initial bordered surface, but the union of these two images in Fig. [6] constitutes 
only the part of the corresponding Riemann surface that is above the new closed geodesies (the 
cyan line); two ends of the green geodesies constitute the double helix approaching the new 
geodesic line but never reaching it, and we always have one copy of the initial surface on one 
side of coils of this helix and the other copy - on the other side. 

3 Algebras of geodesic functions 
3.1 Poisson structure 

One of the most attractive properties of the graph description is a very simple Poisson algebra 
on the space of parameters Z a . Namely we have the following theorem. It was formulated for 
surfaces without marked points in [7] and here we extend it to arbitrary graphs with pending 
vertices. 

Theorem 1. In the coordinates (Z a ) on any fixed spine corresponding to a surface with marked 
points on its boundary components, the Weil-Petersson bracket B WP is given by 



where the sum is taken over all three-valent {i.e., not pending) vertices v and V{, i = 1,2,3 mod 3 ; 
are the labels of the cyclically ordered edges incident on this vertex irrespectively on whether they 
are internal or pending edges of the graph. 

The center of this Poisson algebra is provided by the proposition. 
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Proposition 1. The center of the Poisson algebra (|3.ip ) is generated by elements of the form 
Z a , where the sum is over all edges ofT in a boundary component of F(T) taken with multi- 
plicities. This means, in particular, that each pending edge contributes twice to such sums. 

Proof. The proof is purely technical; for the case of surfaces without marked points on bound- 
ary components it can be found in Appendix B in [5]. When adding marked points, it is 
straightforward to verify that the sums in the assertion of the proposition are central elements. 
In order to prove that no extra central elements appear due to the addition process, it suffices 
to verify that the two changes of the part of a graph shown below, 




do not change the corank of the Poisson relation matrix B(T gtSt s)- ■ 
Example 1. Let us consider the graph in Fig. [8l It has two boundary components and two 

4 4 

corresponding geodesic lines. Their lengths, Y ^ an d Y O^i + 2^i), are the two Casimirs of the 

i=l i=l 

Poisson algebra with the defining relations 

{Y l ,Y i „ 1 } = l mod 4, {Z i ,Y i \ = -{Z i ,Y i - 1 } = l mod4, 
and with all other brackets equal to zero. 

3.2 Classical flip morphisms and invariants 

The Z a -coordinates (which are the logarithms of cross ratios) are called (Thurston) shear coor- 
dinates [221 [T] in the case of punctured Riemann surface (without boundary components). We 
preserve this notation and this term also in the case of windowed surfaces. 

In the case of surfaces with holes, Z a were the coordinates on the Teichmiiller space T g H s , 
which was the 2 s -fold covering of the standard Teichmiiller space ramified over surfaces with 
punctures (when a hole perimeter becomes zero, see [8]). We assume correspondingly Z a to be 
the coordinates of the corresponding spaces T g H s in the bordered surfaces case. 

Assume that there is an enumeration of the edges of T and that edge a has distinct endpoints. 
Given a spine T of S, we may produce another spine r a of £ by contracting and expanding edge a 
of r, the edge labelled Z in Fig. to produce r a as in the figure; the fattening and embedding 
of T a in S is determined from that of T in the natural way. Furthermore, an enumeration of 
the edges of T induces an enumeration of the edges of r a in the natural way, where the vertical 
edge labelled Z in Fig. [7] corresponds to the horizontal edge labelled — Z. We say that T a arises 
from T by a Whitehead move along edge a. We also write r Q/ g = (r^)^, for any two indices a, (3 
of edges, to denote the result of first performing a move along a and then along (3; in particular, 
r aa = r for any index a. 

3.2.1 Whitehead moves on inner edges 

Proposition 2 ([3]). Setting §{Z} = log(e 2 + 1) and adopting the notation of Fig. 0/or shear 
coordinates of nearby edges, the effect of a Whitehead move is as follows: 

W z : (A, B, C, D, Z)^{A + cf>(Z),B - <f>(-Z), C + <j>(Z),D - <f>(-Z), -Z). (3.2) 

In the various cases where the edges are not distinct and identifying an edge with its shear 
coordinate in the obvious notation we have: if A = C, then A' = A + 2(f)(Z); if B = D, then 
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Figure 7. Flip, or Whitehead move on the shear coordinates Z a . The outer edges can be pending, but 
the inner edge with respect to which the morphism is performed cannot be a pending edge. 

B' = B - 2(p(-Z); if A = B (or C = D), then A 1 = A + Z (or C = C + Z); if A = D (or 
B = C), then A' = A + Z (or B' = B + Z ). Any variety of edges among A, B, C , and D can 
be pending edges of the graph. 

We also have two simple but important lemmas establishing the properties of invariance w.r.t. 
the flip morphisms. 

Lemma 1. Transformation (|3.2p preserves the traces of products over paths (|2.5p . 
Lemma 2. Transformation (|3.2p preserves Poisson structure (|3.ip on the shear coordinates. 

That the Poisson algebra for the bordered surfaces case is invariant under the flip transfor- 
mations follows immediately because we flip here inner, not pending, edges of a graph, which 
reduces the situation to the "old" statement for surfaces without windows. 

We also have the statement concerning the polynomiality of geodesic functions. 

Proposition 3. All G 7 constructed by (12. 5p are Laurent polynomials in e Zi and e Yj ' 2 with 
positive integer coefficients, that is, we have the Laurent property, which holds, e.g., in cluster 
algebras [TO]. All these geodesic functions preserve their polynomial structures upon Whitehead 
moves on inner edges, and all of them are hyperbolic elements (G~ > 2), the only exception 
where G 7 = 2 are paths homeomorphic to going around holes of zero length (punctures). 

3.2.2 Whitehead moves on pending edges 

In the case of windowed surfaces, we encounter a new phenomenon as compared with the case 

of surfaces with holes. We can construct morphisms relating any of the Teichmiiller spaces 

m 

and T^ 2 with 5 l = {5{, ... ,5^} and 5 2 = {5\, . . . , 5 2 2 } providing n\ = n2 = n and Yl $1 = 

i=i 

E 5 h that 

is, we explicitly construct morphisms relating any two of algebras corresponding 

i=i 

to windowed surfaces of the same genus, same number of boundary components, and with the 
same total number of windows; the window distribution into the boundary components can be 
however arbitrary. 

This new morphism corresponds in a sense to flipping a pending edge. 

Lemma 3. Transformation in Fig.^is the morphism between the spaces T^ sl and Tgp- These 
morphisms preserve both Poisson structures (|3.ip and the geodesic length functions. In Fig. [9] 
any (or both) of Y -variables can be variables of pending edges (the transformation formula is 
insensitive to it). 
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Figure 8. An example of geodesies whose geodesic functions are in the center of the Poisson algebra 
(dashed lines). Whereas G\ corresponds to the standard geodesic around the hole (no marked points are 
present on the corresponding boundary component), the line that is parallel to a boundary component 
with marked points must experience all possible inversions on its way around the boundary component, 
as is the case for G\\. 



Yi 



Y 2 



Z 



Yi 



-2Z) = Y x 



Y 2 + (j>{2Z) = Y 2 



Figure 9. Flip, or Whitehead move on the shear coordinates when flipping the pending edge Z (indicated 
by bullet). Any (or both) of edges Y\ and Y^ can be pending. 



Proof. Verifying the preservation of Poisson relations (13. ip is simple, whereas for traces over 
paths we have four different cases of path positions in the subgraph in the left side of Fig. [91 
and in each case we have the corresponding path in the right side of this figureH. In each of 
these cases we have the following matrix equalities (each can be verified directly) 

Xy 2 LXzFXzLXy 1 = X Y2 LX Yi , 
X Yl RX z FXzRX Yl = X Yi LX 2 FX 2 RX Yi , 

Xy 2 RXy± = -^y 2 RX ^FX ^ RXy- j 

Xy 2 LXzFXzRXy 2 = Xy R^-z z 
where (in the exponentiated form) 

e^=e Y ^l + e- 2Z Y\ = e Y ^ (l + e 2Z ) , e z = e~ z . ■ 

From the technical standpoint, all these equalities follow from flip transformation (|3.2h upon 
the substitution A = C = Y 2 , B = D = Y\, and Z = 2Z. The above four cases of geodesic 

3 We can think about the flip in Fig. |9]as about "rolling the bowl" (the dot-vertex) from one side to the other; 
the pending edge is then "plumbed" on the left and is protruded from the right side whereas threads of all geodesic 
lines are deformed continuously, see the example in Fig. [22] 



Teichmiiller Theory of Bordered Surfaces 



13 



Yi 



Y 2 




Figure 10. Resolution of the inversion process from Fig.[T] We introduce the new dot-vertex reducing 
the inversion to winding around this vertex (blue part of the path in the graph). Each time a path winds 
around the dot-vertex, we set the inversion matrix F. 



functions are then exactly four possible cases of geodesic arrangement in the (omitted) proof of 
Lemma [TJ 

Using flip morphisms in Fig. [9] and in formula (|3.2p . we may establish a morphism between 
any two algebras corresponding to surfaces of the same genus, same number of boundary compo- 
nents, and same total number of marked points on these components; their distribution into the 
boundary components can be however arbitrary. And it is again a standard tool that if, after a 
series of morphisms, we come to a graph of the same combinatorial type as the initial one (dis- 
regarding marking of edges), we associate a mapping class group operation with this morphism 
therefore passing from the groupoid of morphisms to the group of modular transformations. 

Example 2. The flip morphism w.r.t. the edge Z\ in the pattern in (13. 3|) . 



Y 



Z 2 



-2Zi) y + 0(2Zi) 



z 2 



(3-3) 



where Z\ and Z 2 are the pending edges, generates the (unitary) mapping class group transfor- 
mation 



e Z2 e -Z^ 



e Zl -^e Z2 (l + e 



on the corresponding Teichmiiller space Tjg. This is a particular case of braid transformation 



to be considered in detail in Section 



3.3 New graphical representation 

In the case of usual geodesic functions, there exists a very convenient representation in which 
one can apply classical skein and Poisson relations in classical case or the quantum skein relation 
in the quantum case and ensure the Riedemeister moves when "disentangling" the products of 
geodesic function representing them as linear combinations of multicurve functions. However, 
in our case, it is still obscure what happens when geodesic lines intersect in some way at the 
pending vertex. In fact, we can propose the comprehensive graphical representation in this case 
as well! For this, let us come back to Fig. Q] and resolve now the inversion introducing a new 
dot-vertex at a pending vertex inside the fat graph and assuming that the inversion matrix F 
corresponds to winding around this dot-vertex as shown in Fig. fTTJl 

We now formulate the rules for geodesic algebra that follow from relations (|3.ip and classical 
skein relations. They coincide with the rules in the case of surfaces with holes except the one new 
case depicted in Fig. [TTl Note that all claims below follow from direct and explicit calculations 
involving representations from Section [2} 
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3.3.1 Classical skein relation 

The trace relation tr (AB) + tr (AB^ 1 ) — tr A • tr B = for arbitrary 2x2 matrices A and B 
with unit determinant allows one to "disentangle" any product of geodesic functions, i.e., express 
it uniquely as a finite linear combination of generalized multicurves (see Definition [5] below). 
Introducing the additional factor #G to be the total number of components in a multicurve, we 
can uniformly present the classical skein relation as 




We assume in f|3.4[) that the ends of lines are joint pairwise in the rest of the graph, which 
is the same for all three items in the formula. Of course, we perform there algebraic operations 
with the algebraic quantities - with the (products of) geodesic functions corresponding to the 
respective families of curves. 



3.3.2 Poisson brackets for geodesic functions 

We first mention that two geodesic functions Poisson commute if the underlying geodesies are 
disjointly embedded in the sense of the new graph technique involving dot-vertices. Because of 
the Leibnitz rule for the Poisson bracket, it suffices to consider only "simple" intersections of 
pairs of geodesies with respective geodesic functions Gi and G2 of the form 

Gi = tr 1 • • • XlR l X 1 z L 1 X\ ■ ■ • , (3.5) 
G 2 = tr 2 • • • X 2 B L 2 X 2 Z R 2 X 2 D • • • , (3.6) 

where the superscripts 1 and 2 pertain to operators and traces in two different matrix spaces. 

The positions of edges A, B, C, D, and Z are as in Fig. [71 Dots in (|3.5p . (|3.6p refer to arbitrary 
sequences of matrices R 1 ' 2 , L 1,2 , X z ' 2 , and F 1,2 belonging to the corresponding matrix spaces; G\ 
and G2 must correspond to closed geodesic lines, but we make no assumption on their simplicity 
or graph simplicity, that is, the paths that correspond to G\ and G2 may have self- and mutual 
intersections and, in particular, may pass arbitrarily many times through the edge Z in Fig. UJ 

Direct calculations then give 

{G 1 ,G 2 } = ±(G n -G 1 ), (3.7) 

where Gi corresponds to the geodesic that is obtained by erasing the edge Z and joining together 
the edges "^4" and "£)" as well as "£>" and "C" in a natural way as illustrated in the middle 
diagram in (|3.4p ; G H corresponds to the geodesic that passes over the edge Z twice, so it has the 
form tr • • • XqRzRd XbLzLa • • • as illustrated in the rightmost diagram in (|3.4h . These 
relations were first obtained in [12] in the continuous parametrization (the classical Turaev-Viro 
algebra) . 

Having two curves, 71 and 72, with an arbitrary number of crossings, we now find their 
Poisson bracket using the following rules: 

• We take a sum of products of geodesic functions of non(self)intersecting curves obtained 
when we apply Poisson relation (|3.7p at one intersection point and classical skein rela- 
tion (13. 4p at all the remaining points of intersection; we assume the summation over all 
possible cases. 
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(c) (d) 

Figure 11. An example of two geodesic lines intersecting at the dot-vertex. We present four homotopical 
types of resolving two intersections in this pattern (Cases (a)-(d)). Case (d) contains the loop with only 
the dot-vertex inside. This loop is tiF = 0, so the whole contribution vanishes in this case. The 
(green) factors in brackets pertain to the quantum case in Section [4] indicating the weights with which 
the corresponding (quantum) geodesic multicurves enter the expression for the product GfG§. 

• If, in the course of calculation, we meet an empty (contractible) loop, then we associate the 
factor —2 to such a loop; this assignment, as is known [5], ensures the Riedemeister moves 
on the set of geodesic lines thus making the bracket to depend only on the homotopical 
class of the curve embedding in the surface. 

• If, in the course of calculation, we meet a curve homeomorphic to passing around a dot- 
vertex, then we set tr F = into the correspondence to such curve thus killing the whole 
corresponding multicurve function. 

These simple and explicit rules are an effective tool for calculating the Poisson brackets in 
many important cases below. 

Because the Poisson relations are completely determined by homotopy types of curves in- 
volved, using Lemma [3l we immediately come to the following theorem 

Theorem 2. Poisson algebras of geodesic functions for the bordered Riemann surfaces S S)< 5i 
and ^ 9t s 2 that differ only by distributions of marked points among their boundary components 
are isomorphic; the isomorphism is described by Lemma^ 

It follows from this theorem that we can always collect all the marked points on just one 
boundary component. 
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Gi 



n=3 Mr*" 



T.rrir 



n—4 



-i lip 



G 



2.5 



G 13 'l-J 1 'l-J 1 



G 2 



Figure 12. Generating graphs for ^4„ algebras for n — 3,4, . . . . We indicate character geodesies whose 
geodesic functions G^ enter bases of the corresponding algebras. 



3.4 The A n algebras 

Consider the disc with n marked points on the boundary; examples of the corresponding rep- 
resenting graph T n are depicted in Fig. [12] for n = 3,4,.... We enumerate the n dot-vertices 
clockwise, i,j = 1, . . . , n. We then let Gij with i < j denote the geodesic function corresponding 
to the geodesic line that encircles exactly two dot-vertices with the indices i and j. Examples 
are in the figure: for n = 3, red line corresponds to G12, blue - to G23 and green - to G13. Note 
that in the cluster terminology (see [llj ) these algebras were called the ^4 n _2-algebras. 
Using the skein relation, we can close the Poisson algebra thus obtaining for ^3: 

{Gi2,G 2 3} = G12G23 ~ 2 Ci3 and cycl. permut. (3.8) 

Note that the left-hand side is doubled in this case as compared to Nelson-Regge algebras 
recalled in [5]. In the A3, case this is easily understandable because, say, 

G12 = tv LX2z 2 RX 2Zl = e Zl+Z2 + e Zl ~ Z2 + e~ Zl ~ z \ (3.9) 

and this expression literally coincides with the one for the algebra of geodesies in the case of 
higher genus surfaces with one or two holes (see [4]) but the left-hand side of the relation is now 
doubled (the analogous expression for G12 in [1] was the same as in (|3.9p upon the substitution 
Z\ = X\/2 and Z2 = X2/2, but with the X-variables having the doubled Poisson brackets 
{X2,X\} = 2). In higher-order algebras (starting with n = 4), we meet a more complicate case 
of the fourth-order crossing (as shown in the case n = 4 in Fig. [12]) . Using our rules for Poisson 
brackets, we find that those for these geodesic functions are 

{G13, G24} = 2G12G34 — 2G14G23 (3.10) 

(note that the items in the products in the r.h.s. mutually commute). 

It is also worth mentioning that after this doubling that occurs in the right-hand sides of 
relations (|3.8p and (|3.10j) . we come exactly to algebras appearing in the Frobenius manifold 
approach [6]. 

3.5 The .D^-algebras 

We now consider the case of annulus with n marked points on one of the boundary component 
(see the example in Fig. [8] Here, again, the state of art is to find a convenient (finite) set of 
geodesic functions closed w.r.t. the Poisson brackets^]. In the case of annulus, such a set is given 
by geodesic functions corresponding to geodesies in Fig. [TBI 

We therefore describe a set of geodesic functions by the matrix Gij with i, j = 1, . . . , n where 
the order of indices indicates the direction of encompassing the second boundary component of 
the annulus. 



4 Usually we can say nothing about the uniqueness of such a set for a particular geometry. 
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Figure 13. Typical geodesies corresponding to the geodesic functions constituting a set of generators 
of the D n algebra. We let Gij, i,j = 1, . . . ,n, denote these functions. The order of subscripts is now 
important: it indicates the direction of encompassing the hole (the second boundary component of the 
annulus). The most involved pattern of intersection is on the right part of the figure: the geodesies have 
there eight-fold intersection; in the left part we present also the geodesic function Gu corresponding to 
the geodesic that starts and terminates at the same window. 



Lemma 4. The set of geodesic functions Gij corresponding to geodesies in Fig. [13] is Poisson 
closed. 

The relevant Poisson brackets are too cumbersome and we omit them here because one can 
easily read them from the corresponding quantum algebra in formula (|4.13p below in the limit 
as h — > 0. 



3.6 Braid group relations for windowed surfaces 
3.6.1 Braid group relations on the level of Z- variables 

We have already demonstrated in Example [2] a m.c.g. relation interchanging two pending edges 
of a graph. In a more general case of A n -algebra, we have a graph depicted in Fig. [12] and 
another intertwining relation arises from the three-step flipping process schematically depicted 
in Fig. [Ill 

The graph for the A n algebra has the form in Fig. [12] with Y{, 2 < i < n — 2, being the 
variables of internal edges and Zj, 1 < j < n, being the variables of the pending edges and we 
identify Y\ = Z\ and Y n -\ = Z n to make formulas below uniform. 

We let Ri ; i+i denote the intertwining transformation in Fig. [TUfor 2 < i < n — 2 and in Fig. [9] 
for i = 1 and i = n — 1. For the exponentiated variables, these transformations have the form 

- y- -r c V j- -r <=. 



R 



( e Y *-i ) 




e Y < 




< e^ 1 


> = < 


e z > 








. e Zi+1 . 





-i(l + e 2 ^(l + e^)) 
(l + e 2Z l(1 + e y l)2 y 

l+e 2Z »(l+e y i) 2 



l+e 2Z i(l+e Y i) 



2 < i < n 



D 2Zi+Zi+i+Yi 



-Zi-Yi 



;i + e 2 ^(l + e y »)y 



l + e 2^ (1 + e K i) ) 



(3.11) 



and 





f ezl 1 




Rl,2 < 






{ ^ J 





Z2 (l + e - 

-Zi 



2Zi 



3 Y2 (1 + 



(3.12) 
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li-i Yi Y 



i+1 



Zj z. 



i+1 



Zi z. 



i+1 



Figure 14. Three-step flip transformation of intertwining pending edge variables Zi and Zi+i that 
results in the same combinatorial graph. The rest of the graph denoted by dots remains unchanged. 









R \ 


f J: 1 1 









e Zn (l + e 

e -Z n -! 
0*71- 



-2Z n —-\ 



'l + c 



2Z n - 



(3.13) 



The following lemma is the direct calculation using (|3.1ip . (|3.12p . and (|3.13p . 
Lemma 5. For any n > 3, we /lave i/ie braid group relation 

Ri-l,iRi,i+\Ri,i-l = Ri,i+lRi-l,iRi,i+l, 2 < i < n — \. 



3.6.2 Braid group relations for geodesic functions of ^-algebras 

Here we, following Bondal [2], propose another, simpler way to derive the braid group relations 
using the construction of the groupoid of upper-triangular matrices. It was probably first used 
in [6] to prove the braid group relations in the case of A% algebra. In the case of A n algebras for 
general n, let us construct the upper-triangular matrix A 



A 



( 1 Gi,2 
1 










Gl,3 
^2,3 
1 



G\,n \ 
G*2,n 







G n — 1,71 

1 



(3.14) 



associating the entries Gij with the geodesic functions. Using the skein relation, we can then 
present the action of the braid group element -R^j+i exclusively in terms of the geodesic functions 
from this, fixed, set: 



RiA+\A — A, 



where 



Gj,i+1 

Gi 



Gj,ii 



j > i + 1, 
j < i, 



T h3 ~ GijGi^+i — Gj+ij, j > i + 
Gj,i = Gj^Gi^i+i — Gjj+i, j < i, 
Gii+i = Gii+i. 



A very convenient way to present this transformation is by introducing the special matrices 
Ba+i of the block-diagonal form 



/ 1 



B 



i 

i + 1 



\ 



"i,i+l — 1 
1 



1 / 
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Then, The action of the braid group generator Ri.i+i on A is merely 
Ri s i+1<A = Bi t i + iABf i+l 



(3.15) 



with Bj i+1 the matrix transposed to -Bj^+i. The proof [2] of Lemma [5] in this setting is much 
simpler than in terms of the Teichmiiller space variables; moreover, using this approach, we can 
attack another important issue related to the second braid group relation!. 

We consider the action of the chain of transformations i? n _i in i? n _2, n -i ■ ■ ■ Ri 3R1 lA- Note 

(i—l) 

that, on each step, the item G\ i+ { in the corresponding matrix -Bj^+i is the transformed quantity 

(we assume to coincide with the initial Gij in A). However, it is easy to see that for just 

this chain of transformations, G^.^ = G^j +1 = Gi^+i, and the whole chain of matrices B can 
be then expressed in terms of the initial variables Gij as 



B = B n . 



l,nB n -2,n-l ' ' ' 



-B2,3-Bl,2 



/ Gi, 2 

Gl,3 
Gl,n 

V 1 



-1 








\ 


-1 
/ 



/ 1 


^2,3 


G2,4 


• • • G2, n 


Gl,2 


\ 





1 


£3,4 


■ ■ ■ G3 iU 


Gl,3 










1 


G^ n 


Gi,4 












1 


Gl, n 




\o 










1 


/ 



and the whole action on A gives 



A = BAB T 



and we see that it boils down to a mere permutation of the elements of the initial matrix .4. It 
is easy to see that the nth power of this permutation gives the identical transformation, so we 
obtain the last braid group relation. 

Lemma 6. For any n > 3, we have the second braid group relation 

(Rn-l,nRn-2,n-l ' - - -^2,3-^1,2)™ = Id. 

3.6.3 Braid group relations for geodesic functions of D^-algebras 

It is possible to express readily the action of the braid group on the level of the geodesic functions 
Gij, i,j = 1, . . . , n, interpreted also as entries of the n x n-matrix T> (the elements that are not 
indicated remain invariant): 



Ri i+iD = T>, where < 



Gi-\-i t k 


= Gi t k, 


k / i 


i + 1 


Gi,k = 


Gi t kGi t i+\ — Gi + i^, 


k ^ i 


i + l 


Gk,i+i 


= Gk,i, 


k / i 


i + 1 


« = 


Gk,iGi y i+\ — Gi-,%+1, 


k / i 


i + l 


Gi s i+1 


= Gi t i+i, 






Gi+i t i~ 


-1 = Gi t i, 






Gi t i = 


Gi t iGi.i+i — Gj+i^+i, 






, Gi+l,i 


= Gi+i t i + Gi t i+iG ii — 2Gi i iGi-^i ) i-^i. 







(3.16) 



3 This relation was not presented explicitly in [2], so we consider it here in more details. 
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The first braid group relation follows in this case as well from the three-step process, but it can 
be verified explicitly that the following lemma holds. 

Lemma 7. For any n>3, we have the braid group relation for transformations (|3.16p : 

Ri-l,iRi,i+lRi,i-l'D = Ri,i+lRi-l,iRi,i+l'D , 2 < % < 71 — 1. 

Note that the second braid-group relation (see Lemma E]) is lost in the case of Z) n -algebras. 

To present the braid-group action in the matrix-action (covariant) form (|3.15p note that the 
combinations G^j, G, fc, and GkfcGjj have similar transformation laws in (|3,16|) in the case 
where at least one of the indices j and k is neither i nor i + so we can try to construct globally 
covariantly transformed matrices from linear combinations of the above (coefficients of these 
combinations can be different above and below the diagonal). Note that (since the braid-group 
transformation acts on the A n subgroup of D n in the same way as before), the matrices A (13. 14h 
and A T are transformed as in (|3. 15|) : the analysis shows that we also have two new matrices, 1Z 
and S, with the same transformation law as in (|3.15p : 

{Gjj + Gij — Gi^Gjj j > i, 
— Gij + Gi^Gjj j < i, 
j = i, 

(<S)i,j = Gi^Gjj for all 1 < i,j < n, 
where 1Z is skewsymmetric (1Z T = and S is symmetric (S T = S). 

Lemma 8. Any linear combination w\A + W2A T + plZ + aS with complex ui\, W2, p, and a 
transforms in accordance with formula (|3.15p under the braid-group action. 

We postpone the discussion of modular invariants constructed from these four matrices till 
the discussion of the quantum D n braid-group action in Section 14.5.21 



(3.17) 
(3.18) 



4 Quantum Teichmiiller spaces of windowed surfaces 
4.1 Canonical quantization of the Poisson algebra 

A quantization of a Poisson manifold, which is equivariant under the action of a discrete group D, 
is a family of *-algebras A h depending on a positive real parameter h with T> acting by outer 
automorphisms and having the following properties: 

1. (Flatness.) All algebras are isomorphic (noncanonically) as linear spaces. 

2. (Correspondence.) For h = 0, the algebra is isomorphic as a 2?-module to the *-algebra of 
complex- valued functions on the Poisson manifold. 

3. (Classical Limit.) The Poisson bracket on ^1° given by {ai, 02} = lim ^ ai ^ coincides with 
the Poisson bracket given by the Poisson structure of the manifold. 

Fix a cubic fatgraph T g $ as a spine of S 9 ,<5, and let T h = T h (r g $) be the algebra generated 
by Z^, one generator for each unoriented edge a of T g ^, with relations 

[Z*,Z$]=2mh{Z a ,Zf j } (4.1) 

(cf. (|3.ip ) and the *-structure 

(z*y = zi 
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where Z a and {•, •} denotes the respective coordinate functions and the Poisson bracket on the 
classical Teichmiiller space. Because of (|3.ip . the right-hand side of (|4.1|) is a constant taking only 
five values 0, ±2?rih, and ±4nih depending upon five variants of identifications of endpoints of 
edges labelled a and /?. 

All the standard statements that we have in the case of Teichmiiller spaces of nonwindowed 
Riemann surfaces are transferred to the windowed surface case. 

Lemma 9. The center Z h of the algebra T h is generated by the sums ^ over all edges 

a £ I surrounding a given boundary component, the center has dimension s, and the Poisson 
structure is nondegenerate on the quotient T h /Z h . 

The examples of such boundary-parallel curves are again in Fig. [8l Of course, those are the 
same curves that provide the center of the Poisson algebra. 

A standard Darboux-type theorem for nondegenerate Poisson structures then gives the fol- 
lowing result. 

s 

Corollary 1. There is a basis for T h /Z h given by operators pi, qi, fori = 1, . . . , 3g— 3+s+ ^ 6j 

i=i 

satisfying the standard commutation relations \pi,qj] = 2mHSij. 

Now, define the Hilbert space Ji to be the set of all L 2 functions in the ^-variables and let 
each g-variable act by multiplication and each corresponding p-variable act by differentiation, 
Pi = 2irih-£p. For different choices of diagonalization of non-degenerate Poisson structures, 
these Hilbert spaces are canonically isomorphic. 



4.2 Quantum flip transformations 

The Whitehead move becomes now a morphism of (quantum) algebras. The quantum Whitehead 
move or flip along an edge of T by equation (|3.2|) is described by the (quantum) function [3] 

Trh r e~ ipz 

4>(z) EE <j>\z) = -— / ——————dp, (4.2) 

2 J n smh(7rp) smh(7rnpj 

where the contour f2 goes along the real axis bypassing the origin from above. For each un- 
bounded self-adjoint operator Z h on H, (f) h (Z h ) is a well-defined unbounded self-adjoint operator 
on 7i. 

The function 4> h (Z) satisfies the relations (see [3]) 

4> h {z) - <t>\-z) = z, 

2irih 



(j) h {Z + mh) - (j) h {Z - iirh) 
<p h {Z + in) - <j> h {Z - m) = 



l + e-Z' 
2vri 

1 + e- z / h 



and is meromorphic in the complex plane with the poles at the points {ni(m + nh), m,n € Z + } 
and {— 7rz(m + nh), m,n £ 

The function <j) h {Z) is therefore holomorphic in the strip |ImZ| < 7rmin (1, Re h) — e for any 
e > 0, so we need only its asymptotic behavior as Z € ffi and \Z\ — > oo, for which we have (see, 
e.g., 



\Z)\ ]ZHoo = (Z+\Z\)/2 + 0(l/\Z\). 
We then have the following theorem [31 [13] 
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Theorem 3. The family of algebras T h = T (T 9j g) is a quantization of T g H s for any cubic 
fatgraph spine Tg : s ofEg S, that is, 

• In the limit h~ i— > 0, morphism ()3.2[) using (|4.2p coincides with classical morphism ()3.2p 
with (p(Z) = log(l + e z ). 

• Morphism (|3.2p using (|4.2|) is indeed a morphism of *- algebras. 

• j4 /Zip satisfies W% = I, (|3.2p . and flips satisfy the commutativity relation. 

• F/ips satisfy the pentagon relation. 

• The morphisms T h {T) — > T 1//?l (r) given by i— > ^ commute with morphisms (|3.2p . 
4.3 Geodesic length operators 

We next embed the algebra of geodesic functions (|2.2p into a suitable completion of the con- 
structed algebra T h . For any 7, the geodesic function G 7 can be expressed in terms of shear 
coordinates on T H : 



G-y 



tr P Zl ...z n = J^exp J - mj(j,a)Z a I , (4.3) 



where 01^(7, a) are integers and J is a finite set of indices. 

In general, sets of integers {m,j(j, a)}^ 1 6+3s+2 *' 5 may coincide for different £ J; we 

however distinguish between them as soon as they come from different products of exponentials 
e ±Zi/2 - n traces of matrix products in (|4.3p . 

For any closed path 7 on X Si <5, define the quantum geodesic operator € to be 

= I tr P Zl ...z n ^ E ex P U ^ K-(7,«)^ + 2^(7, a)) I , (4.4) 



where the quantum ordering * • x implies that we vary the classical expression (|4.3p by intro- 
ducing additional integer coefficients Cj(7,a), which must be determined from the conditions 
below. 

That is, we assume that each term in the classical expression ()4.3|) can get multiplicative 
corrections only of the form q n , n € Z, with 

q = e~ inR . 

We often call a quantum geodesic function merely a quantum geodesic because quantum 
objects admit only a functional description. 

We now formulate the defining properties of quantum geodesies. 

1. If closed paths 7 and 7' do not intersect, then the operators G0j and Gy commute. 

2. Naturality. The mapping class group MC(T> g $) (|3.2p acts naturally, i.e., for any {G^}, 
W h G MC(E g>s ), and closed path 7 in a spine T g> s of £g,<5, we have 

W\G h 1 ) = G h w{l) . 

3. Geodesic algebra. The product of two quantum geodesies is a linear combination of quan- 
tum multicurves governed by the (quantum) skein relation below. 

4. Orientation invariance. Quantum traces of direct and inverse geodesic operators coincide. 
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5. Exponents of geodesies. A quantum geodesic corresponding to the n-fold concatena- 
tion of 7 is expressed via exactly as in the classical case, namely, 

G^ = 2T n (G*/2), (4.5) 

where T n (x) are Chebyshev's polynomials. 

6. Hermiticity. A quantum geodesic is a Hermitian operator having by definition a real 
spectrum. 

We shall let the standard normal ordering symbol :e ai e a2 ■ ■ ■ e an : denote the Weyl ordering 

gQ-lH ha n j g 

. e a le a 2 ... e ^ : = i + ( ai + ... + an ) + ( ai + ... + an )( ai + . . . + an ) + . . . 

for any set of exponents with Oj ^ —aj for % 7^ j, In particular, the Weyl ordering implies total 
symmetrization in the subscripts. 

We have [3] the proposition, which can be extended to the case of windowed surfaces assuming 
the modification of the "old" notion of graph simple geodesies. 

Definition 1. For a spine T g> g, we call a geodesic graph simple if it does not pass twice through 
any of inner edges of the graph and has at most one inversion at any of pending edges. 

Proposition 4. For any graph simple geodesic 7 with respect to any spine T, the coefficients 
Cj(7, a) in (|4.4p are identically zero, i.e., the quantum ordering is the Weyl ordering. 

Proof. Let us again denote by Y^, i = 1, . . . , 6g — 6 + 3s + #5, the quantum shear coordinates 
of inner edges and by Zj, j = 1, . . . , #5 the quantum shear coordinates of pending edges. But 

the latter always come in the combination X z nFX z n = X 2Z n, so, considering term- by-term the 

i j j 

trace of the matrix product for a quantum graph simple geodesic, we find that we can expand 
it in Laurent monomials in e Y ^ 2 and e z i . It is easy to see that each term e Y ^ I 2 and e z i comes 
either in power +1, or —1 in the corresponding monomial and there are no equivalent monomials 
in the sum. This means that in order to have a Hermitian operator, we must apply the Weyl 
ordering with no additional g-factors (by the correspondence principle, each such factor must 
be again a Laurent monomial in q standing by the corresponding term, which breaks the self- 
adjointness unless all such monomials are unity). Since quantum Whitehead moves must preserve 
the property of being Hermitian, if a graph-simple geodesic transforms to another graph-simple 
geodesic, then a Weyl-ordered expression transforms to a Weyl-ordered expression, and only 
these expressions are self-adjoint. ■ 

Example 3. For the A 3 algebra graph in Fig. \12\ we have exactly three graph simple geodesies 
with the corresponding geodesic functions G\ 2 i G 23 , and G^ 3 given by formulas (|3.9p (which are 
written already in the Weyl-ordered form), and if we consider, for instance, the product 

G^G^ = q- 1 G^ 32 + qG^ l , (4.6) 

where Gi 232 and G^ 3 correspond to respective cases (a) and (b) of resolving crossing of the 
geodesies 723 and 712 near the dot-vertex 2 in Fig. [Til Note that G^ 3 is also Weyl-ordered, 
G 13 = e z ' A+Zl + e Z3 ~ Zl + e ~ Zi ~ Zl whereas 

^1232 = e Zl+2Z2+Zi + e Zl+2Z2 ~ Zi + e Z ^~ 2Z 2-Z- A + e -Z!-2Z 2 -Z 3 + ^2 + ^-2yZ x -Z-, 



apparently is not Weyl-ordered. Product of the same operators as in (|4.6p but taken in opposite 
order gives G\ 2 G\ 3 = qGi 232 + Q~ l Gi 3 , so, introducing the g 2 -commutator [A, B\ q 2 = qAB — 
q~ 1 BA and £ = q 2 — q~ 2 , we have the quantum ^-algebra: 

trih i^ih 1 _ c/^<h \r~<ti s~*fi 1 — tr^ r~*ti 1 — tr^ 

1^23 > (j \2\q 2 — 4^13' 1^13' Lr 23lq 2 ~ £^12 > 1^12 > ^13Jg 2 ~ ^ ( - T 23' 
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4.4 Quantum skein relations 

We now formulate the general rules that allow one to disentangle the product of any two quantum 
geodesies. 

Let Gi and G\ be two quantum geodesic operators corresponding to geodesies 71 and 72 
where all the inversion relations are resolved using the dot-vertex construction (see Fig. [10]) . 
Then 

• We must apply the quantum skein relation^ 




simultaneously at all intersection points. 

• After the application of the quantum skein relation we can obtain empty (contractible) 
loops; we assign the factor — q — q~ l to each such loop and this suffices to ensure the 
quantum Riedemeister moves. 

• We can also obtain loops that are homeomorphic to going around a dot-vertex; as in the 
classical case, we claim the corresponding geodesic functions to vanish, so we disregard all 
such cases of geodesic laminations in the quantum case as well. 

The main lemma is in order. 

Lemma 10 ([3 [21]). There exists a unique quantum ordering f|4.4|) . which is generated 

by the quantum geodesic algebra ()4.7[) and is consistent with the quantum mapping class groupoid 
transformations (|3.2p . i.e., so that the quantum geodesic algebra is invariant under the action 
of the quantum mapping class groupoid. 



4.5 Quantum braid group relation 
4.5.1 Quantum A^-algebra 

We now consider the quantum geodesic functions associated with paths in the A n -algebra pattern 
in Fig. [12J From the quantum skein relation, it is easy to obtain quantum transformations for 
the quantum geodesic functions Gfj. We introduce the .A^-matrix 





( q 



rift. 
^1,2 

q 


rift 

^1,3 • ' 
rift, 

U 2,3 • ' 


rift 
( - T 2,n 










q 


rift 

"n— l,n 



\ 



(4.8) 



associating the Hermitian operators Gfj with the quantum geodesic functions. Using the skein 



relation, we can then present the action of the braid group element i+1 exclusively in terms 



6 Here the order of crossing lines corresponding to G\ and G\ depends on which quantum geodesic occupies 
the first place in the product; the rest of the graph remains unchanged for all items in (14.71) . 
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of the geodesic functions from this, fixed set: Rf i+1 A h = A , where 



-ih 



j >i + l, 

3 < h 

Gfj = qGfjGf i+1 — q 2 Gf +1 j = q 1 G iji+1 Gf t j — q 2 Gf +lj -, j > i + 1, 



lGj^Gf i+1 - q z G^ i+l — q L G? i+1 G l j i - q z Gj i+1 , j < i, 



-Inh 



2r<h 



(4: 



We can again present this transformation via the special matrices B^ i+1 of the block-diagonal 
form 



/ 1 



i + l 



\ 



-ink 



G 
1 



t.t+i 



1 / 

Then, the action of the quantum braid group generator Rf i+1 on A h can be expressed as the 
matrix product (taking into account the noncommutativity of quantum matrix entries) 



R^ i+1 A h - B i i+1 A h (B^ i+1 ) 



(4.10) 



with (-Bfj+i) the matrix Hermitian conjugate to Bf i+l (its nontrivial 2 x 2-block has the form 

qG?^ i 

q 2 



Using the same technique as above, it is then straightforward to prove the 
following lemma. 

Lemma 11. For any n > 3, we have the quantum braid group relations 

2<i<n-l, (4.11) 



r>h r>h T>li r>h nfl r>h 

rL i-l,i- fL i,i+l- rt i— l,i ~ rl i,i+l rl i-l,i- n 'i,i+l, 



(R 



n—l,n- n n—2,n 



r>h r>h \n 
-1 ' " " - ft 2,3- ft l,2j 



Id. 



(4.12) 



4.5.2 Quantum .D^-algebra 

We now quantize the Poisson algebra of geodesic functions Gij corresponding to paths as shown 
in Fig. [13j We have there eight possible variants of nontrivial intersections shown in Fig. [T5l 



The corresponding quantum permutation relations read3 (q 



-iirh 



Case (a) 



Case (b) 



[Gfj, Gm] 



S I ^kj^li ~ ^jk^il + ^jl^ik ~ ^Ij^ki 

+ (q + q l )(Gf l GjjG k i k — GjyGjjG^M ; 

Q^jl^ij ~ Q u ij^jl —si ZLr il — — u u\S*jj) I 

+ e(? + q~ X )G^ + (q- q- l )G%G%- 



7 Deriving these relations requires a tedious combinatorial analysis based on quantum skein relations formulated 
in Section l4~4l 
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3 I 




j I 








Figure 15. Eight cases of nontrivial intersections of geodesies from the set GV,-, i, j = 1, . . . ,n in the 
case of the D n -algebra. 



Case (c) [G%, Gji] — ^\Gj k G^ — G^Gf^j ; 



Case (d) qG%G% - q~ l G h kl G% = £G&; (4.13) 



Case (e) [G% G»] = e((G^ 2 - {G%)^j 



Case (f ) [Gji , G*\ = £ (G%Gf t - GlG% 

Case (g) qGjjG^j — q l G k jGjj = £,G kk , qG^ k Gjj — q 1 GjjG^ k = £,G kk ; 
Case (h) [Gl G kk ] = (q - q~ l ) (G% - G ki ) . 

Although these relations not only contain triple terms in the r.h.s. but also noncommuting 
terms (this is the price for closing the algebra), they nevertheless establish the lexicographic 
ordering on the corresponding set of quantum variables {Gfj}. 

Lemma 12. Permutation relations treated as an abstract algebra postulated by (|4.13j) satisfy 
the commutation Jacobi identities. 



The proof is tedious but straightforward calculations. Note that algebra (|4.13p is self- 
consistent even without relation to geometry of modular spaces; the similar phenomenon was 
already observed in the case of ^4 n -algebras. As regarding the classification of cluster algebras 
in [11], we produced the corresponding algebras for disc and annulus with arbitrary number 
of marked points (in our approach, a punctured disc is just an annulus with one hole of zero 
perimeter; using the isomorphism in Theorem[2]we can move all marked points to one boundary 
component). We do not know however as yet an example of such an algebraically closed set for 
a disc with two punctures (holes); this case deserves a separate investigation. 

We now provide the quantum version of the braid group transformations (|3.16p . They are 
similar to (|4.9p . the only actual distinction is in the corresponding 2 x 2-block of the matrix T>. 
It is however that distinction that prevents us from writing these transformations in a form 
similar to (|4.10p . 
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We therefore have for the quantum braid group transformation for the D n -algebra 



^i+l,k 

G~ h 
j,k 

^k,i+l 

G k,i 



„r^h fih „2/~<h „— i/~<h /~ih „—2fih 

Q ijr i.i+i Lr k.i ~ Q u k,i+i — Q ^k^i^+i ~ 1 u k,i+i^ 



i^+l^k,! 1 ^jfc.i+l ~ 1 ^k,' 
(-in rih, rih rih 

J,i — < l Lr i,i+l Lr i,i ~ 1 — 1 b i,i Lr i,i+l ~~ 1 ^i+l.i+l' 



A; ^ i,i + 1, 
7^ z,i + 1, 
A; / z,i + 1, 
7^ z,i + 1, 



(4.14) 



Lemma 13. For any n > 3, we have the quantum braid group relations (|4.1ip for transforma- 
tions (|4.14|) of quantum operators subject to quantum algebra (|4.13[) . 

Again, it the second identity (|4.12|) is lost in the case of D n algebras. 



4.5.3 Matrix representation for .D^-algebra and invariants 

We now construct the quantum analogues of (|3.17p and (|3.18p . 

Lemma 14. The following four matrices (with operatorial entries), together with all their linear 
combinations, transform in accordance with the quantum braid-group action (|4.10p ,- A h 
(A h )\ K h , andS h , where 



(TZ r 



Gl i + q 2 Gl J -qGlG h 



"J 



J-' 




q ^Gl i + q ^G^l 



3 > h 
3 < h 
3 = i, 



-TV 



(4.15) 



= for all l<i,j<n, 

Remark 1. We took particular forrr|§ (|4.15p of the matrix TZ h because, in the case n = 2, the 
combination 



l,l ( - T 2,2 



Q ^2,1 — Lr 2,2 Lr l,l Q u l: 



is a central element of the (quantum) algebra D 2 ; the other central element is 

Gi, 2 Gji - q\Gl 2 ) 2 - o" 2 (Gt) 2 = G\ x Gl 2 - q-\G\ % ? - q\G\ x f . 

Also, exactly with such a combination, diagonal elements remain zeros acquiring no o-corrections. 
And, eventually, only this matrix possesses the (quantum) cyclic symmetry below. 

A cyclic permutation of indices P : i \— ► i + 1 mod n; j t— ► j + 1 mod n destroys the structure 
of the matrix A. and results in the following transformations for TZ h and S h : 



P : TV 



( o 



V -q~ 2 
/o i 



i 



\ 

i 



/o 



P : S r 



\ 1 



1 



TV' 



V 



1 



1 

1 \ 



/ 



3 Recall that we can "play" with coefficients adding matrices A n and (.4 
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These transformations together with (|4.10p generate a full modular group. This means that, 
at least in the classical case, det 1Z is the mapping-class group invariant and lies therefore in 
the center of the Poisson algebra. Same is true for S, but det S = whereas det 7Z is nonzero 
for even n = 2m (and vanishes for odd n): denoting Qij := (72.)ij for i < j, we have that 
det 1Z = Pf (TV) 2 , where the Pfaffian Pf(7£) is given by the Grassmann- variable integral 



and is described by all possible (signed) pairings in the set 6>i#2 • • • #2m-i#2m where (6i6j) = Qij 
for i < j. 

For example, for m = 2, we have 

h = Ql,2Q3,4 + Ql,4Q2,3 — Ql,3Q2,4 

(recall that Qij = Gij+Gj^ — Gi^Gjj in the classical case). In the quantum case, these elements 
obviously get ^-corrections to be calculated explicitly since the notion of a quantum determinant 
is ambiguous; we hope to return to this problem elsewhere. 

Also, this construction provides just one central element of the algebra Di m \ finding other 
central elements in a regular way (similar to the one in [2] where all the central elements of the 
-A^-algebra were generated by det(A -1 .A + A.4 T )) needs further investigation. 

5 Multicurves and the double of windowed surface 
5.1 Multicurves for bordered surfaces 

It is a standard notation that a multicurve (lamination) is a collection of non(self)intersecting 
curves. Apparently, in a new formulation with dot-vertices, this definition can be literally 
transferred to the case of surfaces with marked points on boundary components. 

Definition 2. Consider the homotopy class of a finite collection C e = {71, . . . , 7 n } of disjointly 
embedded (unoriented) simple curves 7, in a topological windowed surface ^ g ,s- These curves 
are either closed or terminate at windows. We impose the only restriction that we have an even 
(can be zero) number of endpoints of these curves at each window. An even-based generalized 
multicurve (eGMC) C e in T, g ^ is a multiset based on C; we then have Sj > 1 parallel copies of 
components of C, or in other words, positive integral weights Sj on each component of C, where 
Si is the multiplicity of % in C e . Further, given a hyperbolic structure on £ ffi ,5, we associate to C e 
the product 



of geodesic operators (j2.2j) of all geodesies constituting a GMC; these operators Poisson commute 
in the classical case since the components of C are disjoint. 

We therefore extend the standard definition by introducing nonclosed curves that terminate 
at the boundary components. The only restriction we impose is the evenness of the number 
of curves terminating at each connected part of a boundary component (passing through an 
inversion geodesic line). Now in order to obtain the multicurve, we connect pairwise the ends 
of these curves (below the inversion, or bounding, geodesic line) as shown in Fig. 1161 obtaining 
therefore the collection of topologically closed curves. 

Each curve that is homeomorphic to the path encircling a single dot-vertex (see, e.g., case (d) 
in Fig. 1 11 p has zero geodesic function G = trF, and we therefore exclude all such multicurves. 
At the same time, as before, a eGMC containing a contractible component (of length zero) is 
minus twice the eGMC with this curve removed. 




E 

dBi ■ ■ ■ d8 2m e 1 ^<^ 2 



n 
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Figure 16. Closing the multicurve on a boundary component (shown by blue curves). We let the bullet 
denote the dot-vertex. Lower (yellow) part is fictitious and is not present in the original geometry of 
surfaces with boundaries. 

5.2 The double of the Riemann surface 

In order for our procedure to have a geometrical sense, we now construct the double of the 
Riemann surface and transfer to this double the structure of multicurves. Note that it is often 
easier to prove properties of the geodesic function using the structure of the doubled Riemann 
surface without boundary components. 

We first describe the doubled surface itself. Nonintersecting curves must remain noninter- 
secting in the picture of double when we replace inversions merely by (doubled) variables 2Zj 
on the corresponding edges. To attain this, let us consider the example of a geodesic line in 
Fig. [m If we clone this surface, drift apart the copy of it and then connect all the pending 
ends with their twins preserving the orientation (as shown in Fig. fT8|) . then we obtain the new 
Riemann surface without pending edges (and, correspondingly, without windows). Having the 
initial Riemann surface H g g where S is from (|2.3p . s the number of boundary components, 
and Sj the number of marked points on the jth component (can be zero), we obtain that upon 
the joining as shown in Fig. [TBI each boundary component with even Sj (including the zero one) 
generates two boundary components without marked points (holes) in the doubled Riemann 
surface, whereas each boundary component with odd Sj generates exactly one hole in the new 
doubled Riemann surface. Then, easy calculation using the Euler characteristic formula gives 
the answer for the genus g and number of holes s of the doubled Riemann surface £<?,s : 



Note that while flip morphisms on inner edges of the initial graph pertain to flipping simulta- 
neously two (disjoint) copies of this edge in the double graph, flip morphisms on pending edges 
(see Fig. [9j) change the topological structure of the double graph (the genus and the number of 
holes of the doubled Riemann surface may then change; only the total Poisson dimension must 
remain invariant). 

As in the case of dual Teichmiiller spaces in [HJ, the coordinates Z a describe the linear subspace 
of the Teichmiiller space TjF~ : which is again a 2*-fold covering of the Teichmiiller space ramified 
over punctured surfaces (with old and new boundary components). This subspace comprises 
surfaces admitting the involution interchanging two halves of the doubled Riemann surface. 

It seems that no closed geodesic line corresponds to an initial closed geodesic with odd number 
of inversions (the corresponding line on the doubled Riemann surface is nonclosed). To under- 
stand how to overcome this trouble and to formulate the general procedure of constructing the 
"doubled" multicurve, let us consider the example in Fig. \W\ There, we have the multicurve 
with two components, geodesic curves I and II, on the windowed Riemann surface. Curve I 




s = 2s- #{odd Sj}. 



(5.1) 
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Figure 17. An example of geodesic curve with even number of inversions. 



experiences two inversions while curve II has one inversion. To produce a multicurve on the 
doubled surface, we first merely double the pattern (first stage in Fig. [T9|) and then connect 
each pending edge of the original graph with its copy on the clone of this graph in a natural 
way, that is, preserving the orientation and without introducing intersections of threads of the 
multicurve (second stage in Fig. [T9]) . 

We note, first, that the pattern in the doubled surface thus obtained is always a multicurve 
(since we did not introduce any new intersection). We now turn to its content. 

It is easy to see that a geodesic line I with the geodesic function G\ and the length £j on 
the original surface that has even number of inversions produces two disjoint nonselfintersecting 
(albeit not parallel) geodesic lines Ii and I2 on the doubled Riemann surface, and these two new 
lines has the respective geodesic functions G\ x and G\ 2 and lengths £^ and £\ 2 such that 

G h = G h = G h £ h = £ l2 = h. (5.2) 

On the contrary, for a geodesic line II with the geodesic function G\\ and the length £\\ 
that has an odd number of inversions, this geodesies produces one nonselfintersecting geodesic 
line IIi ; 2 on the doubled Riemann surface with the geodesic function Gn 1 2 and length 2 i an d 
it is easy to see that they satisfy the relations 

Gu ia = (Gh) 2 - 2, £ lll2 = 2£ lh (5.3) 

that is, we have then a single geodesic of doubled length. 

From (|5,2p and (|5.3p it follows that we must take as the characteristic of a multicurve, or 
a rational lamination, the sum of lengths of the constituent geodesies. We must be able to 
construct this sum from the geodesic function or from a quantum geodesic function ensuring the 
positiveness property: the length must be nonnegative function in the classical case or positive- 
definite operator in the quantum case. This is ensured by the following construction proposed 
in[5]. 

Definition 3. The proper length p. 1.(7) = £~ f of a closed curve 7 in the classical or quantum 
case is constructed from the quantum geodesic operator G^ as 

p.l.( 7 ) =L t = 2 lim -log2T n (G^/2), 
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Figure 18. The double representation for the Ricmann surface depicted in Fig. [17] Pending propagators 
on boundary components are connected, each carrying the doubled shear coordinate 2Zi, in a crosswise 
manner. We then obtain an oriented surface whose genus is given by (|5.ip . To each nonself intersecting 
geodesic line that has an even number of inversions in the original picture, we may set into the correspon- 
dence a unique (nonsclfintersccting) line on the doubled Ricmann surface (dashed line for the geodesic 
line in Fig. fT7)) . 

where we take the principal branch of the logarithm and T n are Chebyshev's polynomials 
(cf. (|4.5j) ). Since T n (cosh|) = cosh^, it follows that p. 1.(7) is the hyperbolic length of 7 
in the Poincare metric in the classical case. 

In the operatorial case, we can determine p. 1.(7) explicitly in terms of the spectral expansion of 
the operator G5, which is known exactly [13] • Namely, the basis of eigenfunctions of G0j is labeled 
by the positive number S whereas the eigenvalue corresponding to the eigenfunction \as) has the 
form e s l 2 + e~ s / 2 and the both limits S — > +00 and 5 — > are singular in the functional sense. 
Eigenfunctions \as) constitute an orthogonal and complete basis, so we can define the proper 
length operator to be the one with the same eigenfunctions as and with positive eigenvalues 5/2. 
The operator p. 1.(7) is then a well-defined operator on any compactum in the function space H. 

The proper length of a QMC or GMC C, which we denote as p.l.(C) = £g, is the sum of 
the proper lengths of the constituent geodesic length operators (or the sum of geodesic lengths 
calculated in the Poincare metric in the classical case) weighted by the number of appearances 
in the multiset. 

Lemma 15. Given a general even-based multicurve C e , its proper length 



satisfies the general relation: denoting by Cip the GMC generated on the doubled Riemann 
surface, we have 




irrespectively on the content of this multicurve: how many geodesic curves constitute the corre- 
sponding multicurve, what is the number of inversions experienced by individual curves, etc. 



n 




(5.4) 



i=l 
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Figure 19. A general procedure of doubling a multicurve. We first produce two copies of the original 
pattern, then sew together pairwisc (with orientation preserving and with parallel transport of the multi- 
curve threads) copies of dead ends forgetting about dot-vertices. In the doubled surface, we then obtain 
two copies (Ii and I2) of those original geodesic lines (I) that experienced even number of inversions in 
the original pattern (although these copies arc not parallel, their lengths coincide in the new pattern), 
whereas each original geodesic (II) that experienced odd number of inversions, generates one new geodesic 
(Hi,2 ) of doubled length in the doubled surface. 

We therefore see that the only object for which the Thurston theory can be elaborated are 
sums (|5.4|) of lengths of curves entering a multicurve. Note that we have the convergency 
theorem [5] only for this characteristic. 

But if we want to consider algebraic structures in a consistent way, it is better to do in the 
original pattern (as was done above) because geodesic algebras on a doubled surface do not 
satisfy the standard relations due to the additional involutional symmetry. 

5.3 Elements of Thurston theory of bordered surfaces 

We now extend the notion of train tracks and foliation-shear coordinates to the case of bordered 
surfaces. (In this section, we use the terminology of Thurston's school.) For a comprehensive 
research into the subject, see [2Dj; an introductory description close to the one in the present 
paper is contained in [5]. 
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Figure 20. Freeway from fatgraph. 



5.3.1 Decorated measured foliations and freeways 

We recall material from [18]. If T g ^ C S 9j 5 is a cubic fatgraph spine of S g ,5, then we may 
blow-up each three- valent vertex of T g g into a little trigon as illustrated in Fig. [201 The resul- 
ting object t = Tr s has both a natural branched one-submanifold structure and a fattening, 
and furthermore, components of S 9 ,5\t are either little trigons, or once-holed (once-punctured) 
nullgons, or bigons with one marked point on one of the edges (we therefore forbid monogons). 
Thus, r is not a train track, but it is almost a train track, and is called the freeway associated 
to r 5>< 5. Notice that each edge of Tgj gives rise to a corresponding large branch of r, and each 
vertex gives rise to three small branches. It is easy to see that every measured lamination in 
S 9 < 5 is carried by the freeway r. 

The frontier of a once-holed nullgon component of T,g_s\r is a puncture-parallel curve called 
a collar curve of £ 5i< 5 (curve I in Fig. [8]). We also introduce collar curves that are unions of bigon 
edges constituting a path homeomorphic to a boundary component (curve II in Fig. [BJ). We 
have exactly s nonhomeomorphic types of collar curves. A small branch is contained in exactly 
one collar curve, while a large branch may be contained in either one or two collar curves (and 
each pending branch is contained in exactly one collar curve). 

A measure on a freeway r is a function n € M* B ^ (B(t) is the set of all branches, long and 
short) satisfying the natural switch conditions provided that for each switch u of r, we have 

£ /,(&)= £ 

outgoing incoming 
half —branches b half — branches b 

Such a function satisfying the switch conditions is commonly called a (transverse) measure 
on r, and r itself is said to be recurrent train track if it supports a positive measure \x with 
fj,(b) > for each branch of r. 

For freeways, however, the measure is not necessarily nonnegative (as it is for train tracks). 
Let U(t) denote the vector space of all measures on r. Notice that fj, £ U(t) is uniquely 
determined by its values on the short branches (the switch conditions are equivalent to the 
"coupling equations" 

/i(ai) + fi(pi) = //(e) = //(a 2 ) + //(fo), 

for any long branch e whose closure contains the switches V\ ^ v 2 and dj , bi are the short 
branches incident on Vi for i = 1,2). And, simultaneously, the values on the long branches alone 
also uniquely determine fi. At the same time, the both sets of transverse measures on short 
and long branches are subject to constraints imposed by the positivity condition in the case of 
recurrent train tracks (the triangle inequality in the case of measures on long branches), so the 
problem is to propose new variables describing the transverse measures that are unconstrained. 
We then identify U(r) ~ ]R LB ( T \ where LB(t) denotes the set of long branches of r. 

Recall that, having nonnegative transverse measures on edges of a train track r, we can 
canonically associate to r the corresponding measured foliation JWFo(Yj gi s), which is determined 
up to the equivalence relation that claims that if two foliations differ only by their "collar 
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weights", i.e., the transverse measures of boundary-parallel arcs, then they are equivalent, and 
in each equivalence class we find a unique measured foliation with all collar weights equal to 
zero. 

A nonnegative measure on r then canonically determines a point of MJ-q(F) together with 
a nonnegative "collar weight", and we can define the space MJ r o(F) = A4J r o(F) x W of 
decorated measured foliations. 

We then have the theorem. 

Theorem 4 ([IS]). The space U(t) ~ R LB ( T ) gives global coordinates on MTq, and there is 
a canonical fiber bundle II : MJ-q — > AdJ-Q, where the fiber over a point is the set M s of all 
collar weights on F . 

Remark 2. The natural action of mapping class group MCg is by bundle isomorphisms of II. 

Furthermore, II admits a natural MCJ-invariant section a : MTq — ► M.Tq which is determined 

by the condition of identically vanishing collar weights. The restriction of a to MJ-q C M.J-Q 
gives a piecewise-linear embedding of the piecewise-linear manifold MTq into the linear manifold 
(vector space) MF ~ U(t) « R LB ^ . 

5.3.2 Shear coordinates for measured foliations 

We now give an equivalent parametrization of measured foliations in terms of "Thurston's shear 
coordinates" that are closest analogues of Thurston's shear coordinates Z a on T (X). 

We assign a corresponding signed quantity (positive for right, negative for left) as follows. 
Given a measure [i on the long branches of the freeway r associated to the fat graph spine r C X, 
define the (Thurston's foliation-) shear coordinate of the edge indexed by a to be 

C a = \{(i{A) - (i{B) + M (C) - (i(D)), (5.5) 

in the notation of Fig. [7] for nearby branches. 

If ath edge is a pending edge, we use also formula (|5.5|) but with, say, (jl(C) = (i(D) = 
because the edge a has then just one switch v\ at its upper end. 

From the very definition, ( a is independent of collar weights. 

The shear coordinates Q a are not independent, but the only restrictions imposed are 

2> = o 

for the sums over edges a £ I surrounding any given boundary component; we call these condi- 
tions the face conditions for shear coordinates. Thus, the space of foliation-shear coordinates is 
of dimension LB{r) — s and coincides, in particular, with the dimension of the nondegenerate 
Poisson leaf. For any assignment of shear coordinates, there is a well-defined point of A4J-Q 
realizing them, thereby establishing a homeomorphism between MJ-q and this sub-vector space 
^lb(t)-s q ^lb(t) q £ ghggj. coordinates on the long branches of r. 

Remark 3. Note that for arbitrary integer measures (i(Z a ) on long edges, the foliation-shear 
coordinates (, a on pending (long) edges can be half-integer! They are integers only if we claim 
the evenness condition. In the picture of the double, these coordinates are halves of the corre- 
sponding would-be foliation-shear coordinates Q a on X^ (where the foliation is induced from 
the one in E 9i< 5 in a natural way described in Fig. I19p . That we take halves of them follows from 
the doubling of the corresponding shear coordinate Z a . 

For example, the value of the foliation-shear coordinate on the pending edge for the pattern 
in Fig. [16] is +2; it is half-integer only if we have originally nonclosed curves. 
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n(C) 



(case Cz > 0) 
Figure 21. Foliation-shear coordinates. 
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2max(/j(a),/i(6)) — //(c 



Figure 22. Mapping class group transformation of flipping the pending edges acting on transverse 
measures. 

We now describe the action of the mapping class group (the Whitehead moves) on foliation- 
shear coordinates, that is, derive the analogue of formula (|3.2p for measured foliations, which is 
an elementary calculation using the formulas for splitting. 

Lemma 16. Under the Whitehead move in Fig. [7J the corresponding foliation- shear coordina- 
tes of the edges A, B, C , D, and Z situated as in that figure are transformed according to 
formula ()3.2[) 



with 



M Z : (UXb,(c,(d,(z) 

» (U + <MCz), Cb - <M-Cz), Cc + MCz), Cd - ^h(-Cz), -Cz) 

4>h((z) = ((z + Kz\)/2, 



(5.6) 



(5.7) 



i.e., 4>h(x) = x, for x > 0, and zero otherwise. All other shear coordinates on the graph remain 
unchanged. Formula (|5.6f) holds irrespectively on whether some of edges A, B, C , and D are 
pending edges. 

When flipping a pending edge as shown in Fig. [9l the corresponding foliation- shear coordinates 
of the edges Y\, Y%, and Z undergo the transformation 

M z : (( Yi Xy 2 Xz) >-> (Cyi - <M-2Cz), Cy 2 + <fe(2c», -Cz) 

with 4>h(x) from (|5,7p . The corresponding transformations for the transverse measures are de- 
picted in Fig. \ 
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Remark 4. Comparing expressions for the classical function (p(x) = log(l + e x ) and (|5.7p . one 
finds that the latter is a projective limit of the former: 

4>h(x) = lim — (f)(Xx) = lim — (f) h (\x), 

A^+oo A A— >+oo A 

that is, all three transformations coincide asymptotically in the domain of large absolute values 
(or large eigenvalues for the corresponding operators) of Teichmuller space coordinates {Z a }. 
This property was crucial for constructing quantum Thurston theory of surfaces with holes in [5] . 
Presumably, the corresponding construction can be transferred almost literally to the case of 
bordered surfaces. We do not however consider it in this paper postponing the description of 
quantum Thurston theory of bordered surfaces for future publications. 

6 Conclusion 

We have proposed a (novel) graph description of moduli spaces of Riemann surfaces with windows 
(marked points on the boundary) and demonstrated that it provides all the necessary ingredients, 
including the construction of the double Riemann surface with holes (and without windows). 

The approach in this paper might be close to the one in [11] where mutations pushing marked 
points from one boundary component to another have been considered as well. There, however, 
the corresponding points can be transferred only in pairs, which apparently differs from the 
procedure in this paper. Establishing a correspondence between these two approaches definitely 
deserves further studies. 

Worth mentioning is the open/closed string diagrammatics proposed recently in [15j. Wc 
hope that introducing new variables (on pending edges) may allow attaining a comprehensive 
quantitative description of string theories and, in its possible quantum version, may provide a 
bridge to a string field theory description of the related objects. 

Also note the algebras A n and D n , which play an important role in classification of finite- 
type cluster algebras |llj . Eventually, it seems interesting to generalize the construction in this 
paper to higher Teichmuller spaces introduced by Fock and Goncharov [9]. 
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